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ABSTRACT 

The observed near infrared background excess over light from known galaxies is com- 
monly ascribed to redshifted radiation from early, very massive, PopIII stars. We show 
here that this interpretation must be discarded as it largely overpredicts the number 
of J-dropouts and Lycv emitters in ultra deep field searches. Independently of the de- 
tailed physics of Lya line emission, J-dropouts limit the background excess fraction 
due to PopIII sources to be (at best) < 1/24. As alternative explanations can either 
be rejected (e.g. miniquasars, decaying neutrinos) or appear unlikely (zodiacal light), 
whereas the reality of the excess is supported by the interpretation of the angular 
fluctuations, the origin of this component remains very puzzling. We briefly discuss 
possible hints to solve the problem. 
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1 INTRODUCTION 



Recent measurements of the Near Infrared Background 
(NIRB) (see Hauser & Dwek 2001 for a review) have shown 
an intensity excess with respect to observed light from galax- 
ies in deep field surveys (Madau & Pozzetti 2000; Totani 
et al. 2001). The discrepancy is maximal at 1.4 /im, cor- 
responding to 17 — 48 nW m -2 sr _1 (about 2-5 times the 
known galaxy contribution) (Matsumoto et al. 2005); signi- 
ficative discrepancies are found also at longer wavelengths. 
The large uncertainty on the amplitude of the excess is due 
to the subtraction of the zodiacal light (i.e. sunlight scat- 
tered by the interplanetary dust) contribution. Two differ- 
ent models for the zodiacal light have been so far proposed, 
which predict a higher (Wright 1998), or a lower (Kelsall 
et al. 1998) contribution of this component to he NIRB. It 
has to be noted, though, that the excess cannot be com- 
pletely explained by a combination of galaxy and zodiacal 
light alone, unless current zodiacal light models are seriously 
incorrect. 

Estimates based on theoretical models suggest that the 
most likely origin of such excess is redshifted light from 
early, very massive (M > 100A/q) PopIII stars (Santos, 
Bromm & Kamionkowski 2002; Salvaterra & Ferrara 2003 
[SF]; Fernandez & Komatsu 2005) if these form efficiently 
down to z w 8. The same stars can also account for the ob- 
served small scale angular fluctuations detected in the same 
bands (Magliocchetti, Salvaterra & Ferrara 2003; Kashlin- 
sky et al. 2004; Cooray et al. 2004). In this paper, we review 
this conclusion in the light of the new Spitzer satellite and 



ground-based ultra deep searches for high redshift galaxies. 
Throughout the paper we use the AB magnitude system*. 



2 BASIC EQUATIONS 

The mean specific background intensity, J(^o,zo), at red- 
shift zo observed at frequency vo, produced by a population 
of sources (in our case, the PopIII stellar clusters) charac- 
terized by a comoving emissivity t v {z), can be written as 



J(v , z ) = 



(l + *o) 3 

47T 



tv(z)e 



az 



(1) 



where v — vo(l+z) /(1+zo), dl/dz is the proper line element, 
T e / / is the IGM effective optical depth at vq between redshift 
zo and z (see Sect. 2.2 of SF for a full description of the IGM 
modelling) . 

The source term e„ can be written, in general, as a con- 
volution of the light curve of the stellar cluster l v (t) with the 
source formation rate (per unit comoving volume). Under 
the approximation that the source formation rate is constant 
over the source lifetime n / , and considering the average spe- 
cific luminosity per unit mass, I„, over such timescale, the 
emissivity is given by 



e u (z) 
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dM h 
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* AB magnitudes are defined as AB= —2.5 log(.F„ ) — 48.6, where 
F VQ is the spectral energy density within a given passband in units 
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Hz 
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where n(M^, z) is the comoving number density of halos of 
mass Mh at redshift z given by Press & Schechter (1974). 
The integral represents the collapsed mass per comoving 
volume contained in dark matter halos with mass above 
M rn i n (z). We consider two different threshold masses de- 
pending on the relevant cooling agents: (i) atomic hydrogen 
(i.e. T v i r > 10 4 K); and (ii) molecular hydrogen as com- 
puted by Fuller & Couchman (2000; hereafter FC). The time 
derivative converts the collapse mass per comoving volume 
into a formation rate and the pre- factor /*(fi(,/fi m )T con- 
verts from total mass to stellar mass; /* is the PopIII star 
formation efficiency. 

As for PopIII stars properties, we adopt the light curves 
of massive metal-free stars (Schaerer 2002), including neb- 
ular emission and Lya emission line (Loeb & Rybicki 1999; 
for a detailed discussion, see SF). For simplicity, we con- 
sider a 8— function IMF peaked at 300 Mq, since PopIII 
specific spectra above > IOOMq are essentially identi- 
cal (Bromm, Kudritzki & Loeb 2001). Stellar lifetimes are 
also very similar for very massive stars, so we adopt a single 
value nf = 2 Myr. 

The average flux F uo (z, Mh) from PopIII stellar cluster 
of mass M c i at redshift z is given by 



v ° 4irAv d L {z) 2 



(3) 



where cIl(z) is the luminosity distance, Auo is the instru- 
mental bandwidth, and Vmin^max are the restframe fre- 
quencies corresponding to the observed ones. In Fig. 1 we 
show the expected flux from a PopIII stellar cluster with 
M c i = 10 6 Mq at redshift z s = 10. It is interesting to note 
that such an object could be luminous enough to be detected 
at the NICMOS H 160 = 28 limiting magnitude. 

Eqs. (l)-(2) allow us to compute the diffuse NIRB 
background due to PopIII stars and compare it with the 
data. In order to determine the minimum requirements for 
the sources responsible for the NIR light, we consider the 
maximal contribution of both "normal" deep field galaxies 
(Madau & Pozzetti 2000; Totani et al. 2001) and zodiacal 
light (Wright 1998). We then determine the minimum star 
formation efficiency, /™ , necessary to fit the NIRB spectrum 
(additionally, we assume that PopIII stars do not form below 
Zend ~ 8, as implied by the background spectrum J-band 
peak). We use here the data obtained by DIRBE and the 
recent determination by IRST (Matsumoto et al. 2005) that 
provide important information on the shape of the NIRB 
in the 1.4-4 /im wavelength range, rescaled using the Wright 
(1998) zodiacal light model subtraction. We will refer to this 
solution as the minimal NIRB model, for which /™ = 0.8 
(0.4) for H-cooling (Lb-cooling) halos. Note that this value is 
larger than derived by SF, due to their simplifying assump- 
tion that tj / is of the order of the Hubble time. Finally, for 
the given value of /™, we compute the surface density of 
PopIII stellar clusters required to fit the minimal model in 
the flux range F vo and F VQ + dF VQ as 
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Figure 1. Observed flux from a M c ; = 10 6 Mq PopIII stellar 
cluster at z s = 10. The lines show the total flux (solid), and the 
stellar (dotted), nebular (dashed), and Lyo line (long-dashed) 
contributions. 
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where dV c /dzdfl is the comoving volume element per unit 
redshift per unit solid angle, n c (z, F vo ) is the comoving num- 
ber of objects at redshift z with observed flux in [F VQ , F„ + 
dF UQ ], given by 



nciz ' F "°^dF7^ z ' F ^ Tlf dmJt {Mh ' z) 



(5) 



3 OBSERVATIONAL CONSTRAINTS 
3.1 Spitzer number counts 

Ground-based mid-infrared (MIR) galaxy counts are lim- 
ited by the high thermal background emission and nar- 
row atmospheric transmission window. The InfraRed Array 
Camera (IRAC) instrument on the Spitzer Space Telescope 
provided the first source counts in the wavelength range 
3 < Ao < 10 fj,m. Thanks to the low background available in 
space, IRAC can reach faint limiting fluxes in modest expo- 
sure times. Fazio et al. (2004b) reported the combined data 
from three survey fields of complementary area and depth 
at four IRAC wavelengths, 3.6, 4.5, 5.8, and 8.0 fim, respec- 
tively. The 3.5 and 4.5 /J,m images are affected by source 
confusion at their faint limits; for this reason, Fazio et al. 
(2004b) choose to truncate the source counts at the rela- 
tively bright flux value of 1.7/iJy"t. After star subtraction 
and completeness correction, Fazio et al. (2004b) provided 



t We adopt the 'concordance' model values for the cosmological 
parameters: h = 0.7, C m = 0.3, C A = 0.7, C b = 0.038, <r 8 = 0.9, 
and T = 0.21. 



t In the original paper of Fazio et al. (2004b) data are given in 
instrumental magnitude relative to Vega. Throughout the paper 
we convert these into fluxes for clarity. The flux density of Vega 
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Figure 2. Expected number counts for PopIII sources fitting the 
NIRB excess. The solid (dashed) line refers to H-cooling (H2- 
cooling) halos. Left panel: X = 3.5 fira, right panel: A = 4.5 /jm. 
Data (Fazio et al. 2004b) points are from the Spitzer QSO 
HS1700+6416 field (diamonds) and the Extended Groth Strip 
(triangles). 

extragalactic source counts for the three fields and the four 
IRAC wavelengths. The data for the 3.5 and 4.5 (im images 
are shown in Fig. 2 for the two deepest fields, around the 
QSO HS1700+6416 and in the Extended Groth Strip. Data 
the (widest) field in the Bootes region are not displayed due 
to their higher flux limits. 



Figure 3. Top panel: number of z ~ 10 sources predicted by the 
minimal NIRB model. The selection criterion Jno — Hiao > 1-8 
has been applied and a FOV of 2.6 arcmin 2 has been assumed 
for comparison with NICMOS UDFs. The dotted line shows the 
conservative flux limit for the NICMOS survey of Bouwcns et al. 
(2005). Bottom panel: number of sources with strong Lya line 
emission predicted by the minimal NIRB model. The selection 
criterion J s — NB > 0.3 has been applied and a FOV view of 4 
arcmin 2 has been assumed. The dotted line shows the flux limit 
of the VLT/ISAAC survey of Willis & Courbin (2005). In both 
panels the solid (dashed) line refers to H-cooling (H2-cooling) 
halos. Counts flatten below a certain flux as the threshold mass 
is reached. 



3.2 High redshift galaxy searches 

The detection of galaxies at the highest redshifts is challeng- 
ing because of the redshifting of UV light into the infrared 
and the well-known limitations of current IR instruments. In 
spite of this difficulties, the search for high redshift galax- 
ies, based on the dropout technique, has provided the first 
tentative detection of galaxies at redshifts beyond z ~ 7 
(Bouwens et al. 2004), stimulating a re-analysis of the small 
amount of deep IR data. 

Bouwens et al. (2005) looked at the prevalence of galax- 
ies at z ~ 8 — 12 by applying the dropout technique to the 
wide variety of deep F110W- and F160W-band (hereafter 
Jno and -H160, respectively) fields that have been imaged 
with the Near Infrared Camera and Multi-Object Spectrom- 
eter (NICMOS). The principal data set is constituted by 
two 1.3 arcmin 2 deep NICMOS parallels taken with the Ad- 
vanced Camera for Surveys (ACS) Hubble Ultra Deep Field 
(UDF; each has about 160 orbits data). Complementary 
shallower fields possessing similar Jno + i/ieo imaging are 
analysed as well. The 5a limiting magnitude for the UDFs 
is « 28.6 in Jno and pb 28.5 in Hieo (0".6 aperture). Since 

in the 3.5 /im (4.5 fim) IRAC bandpass is 277.5 (179.5) Jy (Fazio 
et al. 2004a). 



these limiting magnitudes are reached only in the deepest 
regions of the two NICMOS UDFs, we adopted the con- 
servative limiting magnitude Hieo = 28 (Bouwens, private 
communication). The primary selection criterion for high- 2: 
sources is Jno — -f/ieo > 1-8. Using this criterion, Bouwens 
et al. (2005) found eleven sources. Eight of these are ruled 
out as credible z ~ 10 sources, either as a result of detection 
(> 2a) blueward of Jno or because of their colors redward 
of the break (Hieo — K ~ 1.5). The nature of the three re- 
maining sources could not be assessed from the data, but 
this number appears consistent with the expected contami- 
nation from low-redshift interlopers. Hence, Bouwens et al. 
(2005) concluded that the actual number of z ~ 10 sources 
in the NICMOS parallel fields must be three or fewer. 

Willis & Courbin (2005) have presented a deep, narrow 
J-band search (called ZEN) for protogalactic Lya emission 
at redshift z ~ 9, combining a very deep VLT/ISAAC image 
of the Hubble Deep Field South using a narrow-band (here- 
after NB) filter centered at 1.187 /im, with existing deep, 
broad-band images covering optical and NIR wavelengths. 
They searched for z f» 9 Lycv-emitting galaxies displaying a 
significant narrow-band excess relative to the J s band (i.e. 
J s — NB > 0.3) that are undetected at optical wavelengths. 
They identified no sources, matching these requirements, up 
to a NB magnitude limit of 25.2 in a field of ~ 4 arcmin 2 . 
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4 IMPLICATIONS FOR NIRB MODELS 

Having fixed the minimal NIRB model (see Sec. 2), we ask 
what are the predictions of such model in terms of: (i) num- 
ber counts in the 3.5 and 4.5 /tm Spitzer bands, (ii) number 
of J-dropouts that can be detected with NICMOS, and (iii) 
number of strong Lya-emitting galaxies identifications. 

The main results are shown in Figs. 2-3 for the two cases 
of H- or H2-cooling halos. In Fig. 2 we plot the predicted 
differential number counts compared to Spitzer data; Fig. 
3 shows the total number of J-dropouts expected in the 
NICMOS ultra deep observations (top panel) and in the 
ZEN search (bottom panel). 

If the NIRB excess is to be ascribed to PopIII clusters, 
we find that even for the minimal NIRB model a substan- 
tial fraction of the Spitzer number counts is provided by such 
galaxies, particularly at the faintest fluxes. For T vir > 10 4 
K halos, almost all the (3.5 fxrn and 4.5 /jm) IRAC counts 
at the flux limit of the survey should correspond to PopIII 
galaxies, whereas for H2 cooling halos ~ 50% of the ob- 
served sources at 1.7 /xJy should lie at 2 > 8. Even at fluxes 
of 4.4 jiiJy, where confusion noise and completeness prob- 
lems are less important, ~60% (20%) of the detected Spitzer 
sources should be PopIII galaxies in H-cooling (H2 cooling) 
halos. 

The situation becomes even more serious if we consider 
the predictions for the number of 2 ~ 10 PopIII galaxies 
detectable in the NICMOS UDFs. Adopting the same selec- 
tion criterion Jno — Jfi60 > 1.8 as in Bouwens et al. (2005), 
the model predicts that thousands of these sources should 
be detected. More precisely, 1165 (5634) J-dropouts should 
be observed for H-cooling (H2 cooling) halos below the con- 
servative magnitude limit iJieo = 28. These numbers are 
at odd with the mere 3 (tentative) detections reported by 
Bouwens et al. (2005). 

The J-band peak in the NIRB excess has been inter- 
preted (SF; Santos et al. 2002) as due to strong asymmetric 
Lya emission line from PopIII clusters at 2 ~ 9. Such sources 
can be identified in a narrow J-band survey from a signifi- 
cant narrow-band (NB) excess relative to the flux observed 
in a broader J-band. In the bottom panel of Fig. 3 we show 
that, by applying the criterion J s — NB > 0.3, almost 394 
(683) H-cooling (H2-cooling) sources should be identified up 
to a magnitude limit NB — 25.2 in a 4 arcmin 2 field for the 
minimal NIRB model. Instead, the ZEN search implement- 
ing the same selection criterion has yielded no identification. 

Can we reconcile the NIRB-PopIII interpretation with 
the constrain given by the currently available observational 
data ? The Spitzer number counts can be reduced by quench- 
ing the formation of PopIII stars in larger halos as a re- 
sult of the transition to a PopII mode driven by chem- 
ical feedback (Schneider et al. 2002, 2003, 2005). In this 
case the minimal NIRB model is reproduced if PopIII star 
formation occurs in all halos with masses in the range 
M min (FC) < M h < 10 9 M Q with /I™ = 0.8. Chemical feed- 
back would drastically reduce the number of PopIII IRAC 
galaxies without reducing the NIRB intensity, the latter be- 
ing mainly sustained by smaller halos. This solution, unfor- 
tunately, does not solve the problem with the NICMOS fields 
and the narrow J-band search. The reason is well exempli- 
fied in Fig. 1, where we show that at 2 w 10, for the expected 
Lya fluxes, NICMOS would be able to detect stellar masses 



as small as M* = 2 x 10 5 M Q , or M h « 2 x 10 6 M Q . This 
figure is increased to M h « 5 x 1O 7 M for VLT/ISAAC. 
Hence chemical feedback cannot provide a way out. 

We can compute a strong limit on the fraction of NIRB 
light produced by PopIII stars based on the three candidate 
z = 10 galaxies in the UDFs. 

Even if all candidates turn out to be genuine 2 = 10 
PopIII galaxies, the contribution to the NIRB is less than 
0.05(0.04) nW m~ 2 sr" 1 in the J (H) band, i.e. w 1/340 
of the minimum observed NIRB excess. Within these limits, 
the predicted contribution to Spitzer counts at larger wave- 
lengths is completely negligible and no sources should be 
present in the ZEN search. 

Alternatively, if PopIII ionizing photons largely escape 
from star forming environments, the Lya line could be pro- 
duced in the intergalactic medium. The consequent strong 
surface brightness decrease would prevent detection and as- 
sociation with the parent galaxy; however, these Lya pho- 
tons would still contribute to the NIRB. This scenario offers 
a possible explanation for the null ZEN result. However, 
the stellar (as opposed to nebular) emission would make 
2 « 10 PopIII clusters observable as J-dropouts even in 
this case. The three NICMOS candidates (corresponding to 
/I™ = 1.5%) imply then a small NIRB contribution from 
PopIII, amounting only to 1/24 of the minimum observed 
excess. Again, these sources would be below detectability for 
Spitzer at larger wavelengths. 

In conclusion, we find that if the NIRB excess is made 
entirely by high redshift PopIII clusters, these should have 
already been observed in deep galaxy searches in large num- 
bers. The few (possibly none!) detections of such objects 
offers little room for the PopIII explanation of the NIRB 
excess, leaving us with a puzzling question concerning the 
origin of this mysterious background component. 



5 DISCUSSION 

Recent measurements of the diffuse background light in the 
NIR bands have shown an intensity excess with respect to 
the light produced by known galaxies in deep field surveys. 
Theoretical models have suggested that this excess could 
be produced by redshifted light from the first very massive 
(M > 100 M Q ) PopIII stars (Santos et al. 2002; SF) if these 
form efficiently down to 2 ~ 8. The same stars can also 
account for the observed small scale angular fluctuations 
detected in the same bands (Magliocchetti et al. 2003). We 
have reanalysed this conclusion in the light of new Spitzer 
MIR observations, and recent high-2 galaxy searches. 

Our main result is that any model reproducing the min- 
imum value of the excess with emission from PopIII clusters, 
inevitably predicts that a large fraction (in some case all) 
of IRAC (3.5 /xm and 4.5 pm) counts at the faintest fluxes 
is due to galaxies at 2 > 8. Even more seriously, the same 
models overpredict the number of J-dropouts and Lya emit- 
ters identified in deep field searches by a factor of thousands. 
In particular, the interpretation of the background light J- 
band peak as due to redshifted Lya emission from sources 
at 2 w 9 would imply a large number (394/683 depending 
on the main gas cooling agent) of identifications in narrow 
J-band surveys. Instead, no such sources have been detected 
by the ZEN search. 
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The three NICMOS z«10 candidate PopIII galaxies 
limit the contribution of such sources to the NIRB excess 
to be < 1/340. This constraint is somewhat looser if the 
Lya emission originates in the IGM rather than in the star 
forming nebula. In this case the NIRB excess fraction due 
to PopIII clusters is < 1/24. 

As we are facing such a severe problem, it is wise to 
consider possible NIRB data contaminations from spurious 
components. The largest experimental uncertainty is repre- 
sented by zodiacal light subtraction. Although it is in prin- 
ciple possible to explain the excess in the NIRB with zodi- 
acal light, this might imply a revision of the properties of 
dust grains canonically assumed (Dwek, Arendt & Krennich 
2005). On the other hand, the existence of the NIR excess 
above the contribution of deep field galaxies at z < 5 seems 
to be suggested by the analysis of the TeV spectral fea- 
tures of distant blazars revealed by Cherenkov experiments 
(although different interpretations cannot be excluded, see 
Aharonian et al. 2005). Finally, zodiacal light has been ex- 
cluded as a possible source of the observed NIRB angular 
fluctuations (Kashlinsky et al. 2002, 2005). 

The interpretation of the NIRB excess must then be 
different. Are there viable alternatives ? Cooray & Yoshida 
(2004) have speculated that a fraction of the unresolved 
NIRB can be due to high-z miniquasars. However, this hy- 
pothesis can be readily discarded by noting that emission 
from a population of early miniquasars significantly con- 
tributing to the NIRB, largely exceeds the unresolved frac- 
tion of the Soft X-ray Background (Salvaterra, Haardt & 
Ferrara 2005; Madau & Silk 2005). In particular, Salvaterra 
et al. (2005) have shown that the SXRB constrains the con- 
tribution of these sources to be < 10~ 3 nW m 2 sr 1 . High 
redshift supernovae have also be examined by Cooray & 
Yoshida (2004), who exclude any significant NIRB contri- 
bution from such objects. Finally, more exotic sources as 
sterile neutrinos have been firmly excluded as viable expla- 
nations for the excess by Mapelli & Ferrara (2005). 

What are the remaining possibilities ? Local, low sur- 
face brightness galaxies undetected by available instruments 
might reduce the gap between NIRB data and the integrated 
light of deep field galaxies; similarly the light from the outer 
regions of moderate redshift galaxies could be missed (Bern- 
stein, Freedman & Madore 2002). Although the nature and 
properties of such unidentified objects are unknown (note 
that these sources must have J > 29), these local galaxies 
would not explain the excess in the small angular scale fluc- 
tuation observed in the NIRB in the J band (Kashlinsky et 
al. 2002). A more promising alternative would be to locate 
these faint sources at z « 2 — 3, for example if the Lyman 
Break Galaxies luminosity function is extremely steep. 

In any case a considerable amount of study will be nec- 
essary to fully understand the puzzling nature of the NIRB 
excess. 
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